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Abstract 

Surface level instability when tube is injected into vibrating bed of powder, 
which was originally found in experiments, is investigated numerically. We 
find that thicker (thiner) tube makes surface level inside tube higher (lower) 
than surface level outside tube. With fixed acceleration amplitude of vibra- 
tion, surface level inside tube becomes higher as amplitude of vibration in- 
creases, which can be explained by considering the dependence upon strength 
of convective flow. 
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I. INTRODUCTION 



The dynamics of granular material attracted many attentions of physicists [[[]. Among 
them, vibrating bed of powder was studied by Faraday half and a century ago 0, and many 
papers were written about its behavior. This is because vibrating bed can exhibit many 
interesting phenomena like convection ||, surface heaping |BJ, surface fluidization ||, size 
segregation J7[, and turbulence M. Other than those phenomena, surface level instability 
due to injection of tube into vibrating bed of powder was observed in experiments M. The 
main purpose of this paper is to reproduce this phenomenon numerically. 

The organization of this paper is as follows. In Sec. |l| we briefly summarized exper- 
imental findings by Akiyama and Shimomura. Numerical modeling used in this paper is 
explained in Sec. |T| and results will be presented in Sec. |V|. Summary and discussion can 
be found in Sec. |V[ 

II. EXPERIMENTS 

Vibrating bed of powder is a vessel filled with granular matter, typically mono-disperse 
glass beads, and vessel is shaken vertically as strong as gravity acceleration. When acceler- 
ation amplitude of vibration exceeds critical value, which is usually a little bit larger than 
gravity acceleration, the bed exhibits several instabilities, e.g., surface heaping 0, surface 
fluidization ||, and convection ||. Further increase of acceleration amplitude results in 
disappearance of heap, and surface starts to fluctuate violently. In order to measure shear 
friction in this vibrating bed, Akiyama and Shimomura |§ injected thick tube into vibrating 
bed and found vary surprising effect. With tube fixed in space, i.e., tube does not vibrate, 
the surface level inside tube differs from surface level outside tube (Fig. p, even if heaping 
is not observed without tube. Thus, this instability may be different from heaping insta- 
bility observed in vibrating bed. The surface level difference depends upon several physical 
parameters as diameter of tube, acceleration amplitude, and particle diameter. It was rather 
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difficult to understand these dependences because experiment of powder can be affected by 
many other fine differences of conditions, e.g., moisture, temperature, and so on. Thus, 
numerical investigation is much suitable to understand these dependences in detail. 



III. NUMERICAL MODEL 

Although there are many numerical schemes for investigating dynamics of powder, we 
employ here distinct element method (DEM) [HJ to reproduce Akiyama's experiments. In 
DEM, granular particle is models as visco-elastic particle, whose interaction is limited within 
short range. Particularly we employ here non-spherical models introduced by Poschel and 
Buchholtz JTTJ, which is known to reproduce static friction effect better than conventional 



models. In their non-spherical model, each granular particle is modeled as a set of five 
sub-particles (See Fig. ||). Although tangential force is ignored, static friction can be 
considered as interaction among surrounding four small sub-particles. Since relative position 
between center sub-particle and surrounding small sub-particles are considered, rotation of 
non-spherical particle can be considered effectively. Each sub-particle obeys the following 
equation 

d 2 fi 



Y, F lnner + E F outer X 0(4~ | fj ~ ?i |) (1) 



/ y - inner ' / y - ouier " \ "%j 



where i^„ ner and F^ teT are interaction within a non-spherical particle and interaction be- 
tween non-spherical particles respectively. Here subscript % is c or r depending upon whether 
% sub-particle is center (c) sub-particle or surrounding (r) sub-particle, is mass of % sub- 
particle and fi is position vector of i sub-particle. 9 is step function, and i'^ is distance 
between sub-particles which belong to different non-spherical particles, 



2R r (ij = r ,r) 

2R C (i,j = c,c) (2) 
R c + R r = c, r) 



where Ri is radius of i sub-particle. F^^. is defined as 



- 1 inner 



—k(\ Ti — r — j 7 ( — J • riij 



fkj (3) 



with riij = (fj —fi)/ | fj — fi \. k and 7 are elastic and viscosity constant respectively, which 

are related to coefficient of restitution, e = exp I , 7r7 ) . La is distance between 

sub-particles among a non-spherical particle, 

R c + R r = c,r) 

(4) 

V2(R C + R r ) (i,j = r,r) 

This means that five sub-particles have interaction only when their relative positions deviate 
from that shown in Fig.0. F l outer . is defined as 



■'■outer 



-k(\ n-fj I -/.) 7 



dfj df 
dt dt 



riij 



n 



(5) 



In order to simulate vibrating bed of powder, we have to introduce vessel composed of 
sub-particles. Interaction between sub-particles which construct non-spherical particles and 
sub-particles which construct vessels employs the same functional form as F^ J Mter , but i 1 ^ is 
replaced with 

I R c + R v = c,v) 
I Rr + Rv = r,v) 

Motion of sub-particles which construct vessel is not affected by the interaction with sub- 
particles which constructs non-spherical particle, but follows given motion of vessel itself, 
e.g., vibration or static state. In addition to this, 'solid plate' is added to the bottom to 
prevent particles from falling out of the vessel through bottom. This plate causes vertical 
force applied to ith. sub-particle, 

^bottom = k [ z bottom ~ z i\i (7) 

where Zb oUo m and Zi are vertical components of bottom and ith. sub-particle, respectively. 
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IV. THE RESULTS 



In the simulations, we employ following parameters (See Fig.|3|). Radius of sub-particles 
are R c = 3.0, R r = 0.5, and R v = 1.5. We ignore distribution of radius and use identical 
non-spherical particles. Number of non-spherical particles are 200, this means, number of 
sub-particles is 1000. Interaction parameters are taken as k/rrti = 600.0 and 7/771, = 2.0 
independent of the kind of sub-particles. Thus coefficient of restitution e = 0.88. k! is take 
to be 2k. For size of vessel, its diameter D v = 90 and its height H v = 135. Tube whose 
height H c = 10.5 is separated from bottom of vessel by H s = 45. Gravity acceleration is 
taken to be 9.8. 

Control parameters of simulation are the acceleration amplitude of vibration T, angular 
frequency of vibration u, and diameter of tube, D c . Here T is related with u as V = auj 2 / g 
where a is amplitude of vibration. The values of these parameters used in simulations are, 
T = 5.0,6.0, and 7.0, u = 3.0,4.0, and 5.0, D c = 30,24, and 18. Total length of simulation 
is over 40 periods. 

Figure f| shows typical snapshot of simulation. In order to measure surface level differ- 
ences, we define height of surface as follows. For example, surface height inside tube Z{ is 
taken such that number of non-spherical particles above Zj is half as much non-spherical par- 
ticles as aligned along section of tube. (For detail, see Fig. |5|.) The height z outside tube is 
defined in similar way. Then height difference between inside and outside tube Az = zi — z Q . 
In addition to this, we measure flow of non-spherical particles within vessel following ref. 
||, i.e., by dividing vessel into cells as large as non-spherical particle, and flow is defined as 
number of particle transport between cells per period of vibration. 

Figure || shows the dependence of Az upon several parameters, and Figure [7| shows flow 
patterns and strength for several parameter values. 
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V. SUMMARY AND DISCUSSION 



In this paper, we numerically investigated height difference Az between inside and outside 
injected tube. Az increases as diameter D c of injected tube increases except for a few cases 
(See Fig. 0). When accelerated amplitude of vibration T and D c are fixed, Az decreases as 
angular frequency of vibration uj increases. Since convection becomes weak as uj increases 
except for a few cases (See Fig. Az can be regarded as a function of convection indirectly. 
These results are schematically shown in Fig|| Upper raw corresponds to larger uj and lower 
raw to smaller. Left column shows the results for smaller D c and right column shows those 
for larger. Thus, in upper left case, Az takes minimum, and in lower right case, Az takes 
maximum. Also convection which occurs when uj is small enough is drawn in lower raw. 

Possible qualitative explanation of these results are as follows. Granular material inside 
tube is affected by both upward and downward forces. Upward force is due to convection. 
Convection is upward at the center of vessel, thus it push up granular material inside tube. 
Difference between upper raw and lower raw is the difference of strength of convection. 
Lower raw has stronger upward force due to convection, thus surface level inside tube is 
higher in lower raw than in upper raw. 

Downward force is possibly due to difference of density between inside and outside tube. 
Generally, fluidized granular matter has smaller density than that of fixed bed. However, 
granular matter inside tube is hard to flow due to friction with tube, thus density is relatively 
high. This effect pushes down granular matter into tube. This conjecture explains why thiner 
tube has lower level of granular matter inside tube. Thiner tube causes higher friction which 
prevents granular matter from being fluidized, thus heavier granular matter is hard to rise 
by convection. It is schematically illustrated in Fig.[8], where left column has lower level of 
surface inside tube than right column. This conjecture should be confirmed in future. 
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FIGURES 
FIG. 1. Schematics of experiment Q 

FIG. 2. Each non-spherical particle composed of five sub-particles 

FIG. 3. Schematics of numerical simulation 

FIG. 4. Snapshot of simulation. Segments attached to each particle reveals instantaneous 
velocity vectors 

FIG. 5. Definition of height inside tube. Ni{z) is number of non-spherical particles above z. 
The hight Zi inside tube is defined such that N^{z) = m.i/2 where rrii is number of non-spherical 
particles aligned along section of tube (rrii = 6 in this figure). 

FIG. 6. Az. O : co = 3.0, + : u = 4.0, □ : oj = 5.0. (a) T = 5.0 (b) T = 6.0 (c) T = 7.0 

FIG. 7. Flow patterns and its strength (I), (a) T = 5.0, oj = 3.0. Wide tube: / = 1302, 
middle tube:/ = 1103, narrow tube:/ = 1891. (b) V = 6.0, u = 3.0. Wide tube: / = 2821, 
middle tube:/ = 2426, narrow tube:/ = 2546. (c) T = 7.0, u = 3.0. Wide tube: / = 2800, 
middle tube:/ = 2426, narrow tube:/ = 3409. (d) T = 5.0, u = 4.0. Wide tube: / = 446, 
middle tube:/ = 436, narrow tube:/ = 433. (e) T = 6.0, oj = 4.0. Wide tube: / = 686, mid- 
dle tube:/ = 707, narrow tube:/ = 554. (f) T = 7.0, uj = 4.0. Wide tube: / = 931, middle 
tube:/ = 746, narrow tube:/ = 551. (g) T = 5.0, uj = 5.0. Wide tube: / = 368, middle 
tube:/ = 235, narrow tube:/ = 212. (h) T = 6.0, uj = 5.0. Wide tube: / = 311, middle 
tube:/ = 276, narrow tube:/ = 289. (i) V = 7.0, uj = 5.0. Wide tube: / = 257, middle tube:/ = 309, 
narrow tube:/ = 233. 

FIG. 8. Schematics of conclusion 
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